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,-0 . Abstract During the solar journey through galactic space, variations in the physical 

^ ' properties of the surrounding interstellar medium (ISM) modify the heliosphere and 

modulate the flux of galactic cosmic rays (GCR) at the surface of the Earth, with con- 
sequences for the terrestrial record of cosmogenic radionuclides One phenomenon that 
needs studying is the effect on cosmogenic isotope production of changing anomalous 
cosmic ray fluxes at Earth due to variable interstellar ionizations. The possible range 
of interstellar ram pressures and ionization levels in the low density solar environment 
generate dramatically different possible heliosphere configurations, with a wide range 
of particle fluxes of interstellar neutrals, their secondary products, and GCRs arriving 
at Earth. Simple models of the distribution and densities of ISM in the downwind di- 
rection give cloud transition timescales that can be directly compared with cosmogenic 
radionuclide geologic records. Both the interstellar data and cosmogenic radionuclide 
data are consistent with two cloud transitions, within the past 10,000 years and a sec- 
ond one 20,000-30,000 years ago, with large and assumption-dependent uncertainties. 
The geomagnetic timeline derived from cosmic ray fluxes at Earth may require ad- 
justment to account for the disappearance of anomalous cosmic rays when the Sun is 
^ . immersed in ionized gas. 
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1 Introduction 

The Sun has traversed multiple interstellar clouds during it 220 Myr journey around the 
galactic center, including dense neutral clouds, low density partially ionized interstellar 
matter (ISM) such as now surrounds the heliosphere, and hot very tenuous plasma. 
Variations in the cosmic ray fluxes at the surface of the Earth are strongly dependent 
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on the geomagnetic field, the solar magnetic activity cycle, the heliosphere, and the 
physical properties of the circumheliospheric interstellar material. For example, the 
modula tion of the 1 AU cosmic ray flux by solar flare mass ejections has long been 
known |Goslinslll964h . 

The heliosphere acts as a weather vane for the circumheliospheric ISM ( CISM), 
respo nding to the ionization, magnetic pressure, and dynamic ram pressure I Holzerl 
1989). The cosmic ray component at 1 AU varies with the properties of the heliosphere 
modulation region. The interpretation of the geological record of cosmogenic isotopes 
relies on accurate models of the cosmic ray spectra. One factor that is not included in 
the interpretation of the geological record of cosmogenic isotopes is that the cosmic ray 
spectrum incident on the Earth consists of two components that behave differently as 
the Sun travels through space. Galactic cosmic rays dominate at high energies, > 500 
MeV, and are subject to heliospheric modulation as the Sun travels through space. 
However a second cosmic ray component at lower energies is formed inside of the 
heliopause from interstellar neutrals that penetrate and are ionized inside of the helio- 
sphere, forming pickup ions. These are subsequently accelerated to form lower-energy 
anomalous cosmic r ays (ACRs) with a composition derived from neutral interstellar 
atoms in the CISM |Fisk et ffll.lll974r ). The local interstellar cosmic ray spectrum that 
creates the geological radio-isotope record is thus composed of two components that 
vary differently over time and space, the higher energy galactic cosmic rays (GCRs) 
that are modulated by a variable heliosphere, and the ACRs that also depend on the 
density and fractional ionization of the surrounding interstellar cloud. 

In this paper we present the overall picture of the ISM characteristics that result 
from the motion of the Sun and interstellar clouds through space. Observations of in- 
terstellar absorption lines towards nearest stars show that spatial variations in velocity, 
temperature, and ionization of the circumheliospheric ISM create temporal variations 
in the heliosphere boundary conditions. These then cause temporal variations in the 
spectrum and fluxes of cosmic rays at Earth. We also draw possible connections between 
interstellar cloud transitions and the geological radio isotope record. 



2 Contemporary and Paleointerstellar Circumheliospheric ISM 

2.1 Environment, Dynamics, Structure, and Magnetic Field in Nearby ISM 

Interstellar clouds establish the boundary conditions of the heliosphere, and the he- 
liosphere is the primary modulation region for the galactic cosmic ray flux at Earth. 
We therefore briefly review the properties of nearby ISM that sets the future and 
paleointerstellar CISM. 

The Sun is traveling through a region with very low mean densities that extends 
70-200 pc from the Sun (Figure 1).Q Stellar winds and supernovae in star formation 
regions bordering the Local Bubble both sculpted the Local Bubble and contributed to 
the cosmic ray flux at Earth. Hot lo w density plasma is widespread in the Local Bubble 



cavity , T ~ 10 6 K, n < 0.005 cm" 3 jFrisch et aZ.ll2009l . F09). Within ~ 15 pc, a cluster 



1 In Fig. 1, the cumulative dust distributions in and beyond the Local Bubble boundaries 
are evaluat e d from the color excess E(B — V) and astromctric data in the Hipparcos catalog 
(Pcrryman (1997)), after first ignoring stars with variability as shown by the Hipparcos index 
I> 0. The ISM close to Sco-Cen stars, and associated with the 18.5° tilt of Gould's Belt, are 
within the range of |Z| < 50 pc. 
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of local interstellar clouds is flowing past the Sun. Clouds with densities of w 0.3 cm -3 
and « 10 3 cm~ 3 have been identified in this flow. Despite the large difference in scale 
sizes between the heliosphere, with a distance to the upwind heliopause of ~ 150 AU 
(0.0007 pc), and the Local Bubble, the heliosphere traces the solar environment that 
is set by the Local Bubble interior, in particular the interstellar radiation field and 
magnetic field. The Local Bubble environment of the Sun affects the past, present, and 
future heliosphere boundary conditions. 

The heliosphere varies over geologically short timescales due to the velocities of the 
Sun and wispy local interstellar clouds through space. Fig. 1 shows the vector motion 
of the Sun through the interior of the Local Bubble, based on the solar apex motion 
through the local quasi-inertial frame known as the local standard of rest (LSR, we use 
a solar vel ocity of 18.0 ± 0.9 km s" 1 , towards i,b= 47.9° ±3.0°, 23.8° ±2.0°, based on 
results in ISchonrich et cuboid ). 

The heliospheric boundary conditions are set by the CISM. Interstellar neutral gas 
and dust are replenished in the heliosphere over timescales of ~ 30 years due to the 
relative 26 km s _1 motion of the Sun and CISM. Self-consistent photoionization models 
of the surrounding ISM provide a good ma tch to data on the ISM inside and adjacent 
to the heliosphere (|Slavin and Frischlliooj . SO8).0 These models predict densities and 
temperatures for the CISM of n(H°)=0.19 cm~ 3 , n c =0.065 cm~ 3 , n p =0.055 cm -3 , and 
T=6,300 K (from Model 26 in SF08). The fractional ionizations of H and He are 22% 
and 39%, respectively. Constraints on the models include the density and temperature 
of interstellar He inside of the heliosphere, the local ISM towards e CMa (and Sirius), 
pickup ion data giving cloud neutrality, and an interstellar radiation field that includes 
ultraviolet, soft X-ray, and extreme ultraviolet emissions, including contributions from 
a conductive boundary on the local interstellar cloud. Regardless of whether the Local 
Bubble plasma contains hot plasma, the models and data indicate that the CISM is 
warm, partially ionized and low density (see Model 42 in SF08). 

The ISM within ~ 30pc has been shocked, since the gas-phase abundances are 
characteristic of the patte rn expected from the destruction of refractory grains in 50- 
100 km s" 1 shock fronts (|Frisch et a/.lll999h . An association of the local ISM with an 



evolved superbubble shell is thus suggested. The abundances for the refractory elements 
Mg, Si, and Fe are larger by factors of 4-8 than abundances of a cold cloud, and C is 
overabundant. Comparisons of interstellar gas-phase Mg, Fe, and Si abundances with 
solar abundances show tha t the local grains a re olivine silicates, that can be efficiently 
aligned by magnetic fields ( Frisch et qZ.ll201 IT ). 



The properties of the circumheliosphere ISM have varied rapidly over time. LSR 
velocities of the Sun and most nearby interstellar clouds are near 18 km s _1 (or « 18 
pc Myrs -1 ); however the cluster of local interstellar clouds (CLIC, within 30 pc) is 
part of a general flow of ISM past the Sun so that many relative Sun-cloud velocities 
are larger than 14 km s _1 (Table 1). The LSR vector motion of the CLIC is -16.8 
km s _1 arriving from the upwind direction £~ 335°, 6~ —7° (Figure 1). This direc- 
tion is towar ds the cente r of the Loop I (North Polar Spur) supernova remnant (at 
^,6=320°, 5°, lHeiledll998h . and makes an angle of 70° ±35° with the direction of the 



local interstellar magnetic field (jFrisch et al.l 120101 ). These relative Sun-CLIC veloci- 



2 Note there are different velcity models for the Local Interstellar Cloud (LIC, 
lLallement et aZ.lll995l: iFrisch et aDl2002l : iRedfield and Linsky||20u8h : we assume that the ISM 
inside of the heliosphere is part of the LIC seen towards Sirius, 2.7 pc, and e CMa. Recent 
IBEX observations are consistent with the conclusion that the Sun is in the LIC cloud seen 
towards Sirius (see Frisch and McComas, this volume). 
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Fig. 1 Left: The distribution of interstellar dust within 200 pc of the Sun and 50 pc of the 
galactic plane, according to the cumulative amount of ISM traced by color excess E(B-V). 
The solar and CLIC motions through the LSR are given, respectively, by the thick black and 
red arrows (in both figures). The intersection of the SI (cyan) and S2 (gra y) magnetic shell s 
of Loop I with the galactic plane are shown based on spherical shells from Wollcbcn (2007). 
The E(B - V) contour levels of 0.08, 0.11, 0.14, 0.16 mag correspond to N( H°+H 2 ) column 
densities of 4.46e2 6.56e20, 7.87e20, and 9.18e20 cm" 2 , for jV( H°+H 2 )/E(B - V)=5.8e21 
atoms cm -2 mag llBohlin et al\ [I1978T 0. The large blue circles show the three subgroups of 
the Sco-Cen Association; the arc centered near £~ 260° shows the approximate nearside of the 
Gum Nebula. Right: The locations of the nearest stars in Table 1 are plotted in the galactic 
plane , together with the interstellar magnetic field direction (gray arrow, from Frisc hei al\ 
2010), and interstellar clouds within 5 pc of the Sun and 75° of the galactic plane. The clouds 
that are plotted include the LIC and the Blue c louds (green , blue, from RL08), and the G-cloud 
and Apex clouds (lavender and orange, from Frisch 2003). The clouds are labeled, and LSR 
velocities are shown. 

ties and the configuration of local ISM (Figure 1) indicate that the Sun has recently 
emerged from the hot Local Bubble plasma and entered the CLIC. 

Cloud velocity is a key interstellar variable affecting the heliosphere boundary con- 
ditions, because the interstellar ram pressure varies as the square of the Sun-cloud 
velocity. The CLIC is a decelerating flow. From upwind to downwind, cloud velocities 
relative to the Sun (Vkr:) are -28.4 km s -1 t owards 36 Oph, 26.3 km s _1 in the inner 
heliosphere (according to iMobius et qZ.|[20o3 ) ■ and 23.4 km s _1 towards x 1 Ori (Table 
1). For otherwise similar clouds, these velocity differences lead to a 50% difference in 
the ram pressures of the ISM on the heliosphere over timescales of less than « 4000 
years. Over longer timescales of 10 5 years, the heliocentric cloud velocities and cloud 
lengths in Table 1 suggest variations of a factor of ~ 4 in interstellar ram pressures. 

The distribution of the low column density local clouds, A r (H°)< 10 18 ' 5 cm -2 , 
is patchy. The volume density of neutral gas in nearby clouds is known only for the 
CISM. If all local clouds have this same density, n(H°)=0.2 cm -3 , then ~ 35% of the 
sightlines to stars within 10 pc are filled with warm low density ISM. The clouds have a 
mean thickness 0.9 ±0.3 pc, and the Sun would cross these clouds with a mean crossing 
time of ~ 47, 000 years. 

Magnetic fields permeate the region of the CLIC, and create asymmetries in the 
heliosphere configuration. The Sun appears to be located in or near the rim of the 



"SI" magnetic shell associated with Loop I (e.g. IWolleb"er]|2007l : lFriscbll201(f l. Magnet- 
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ically aligned interstellar dust grains near the Sun create a birefringent medium with 
lower optical opacities for directions parallel to the interstellar magnetic field (ISMF), 
and therefore polarize starlight. The direction of ISMF over the nearest 40 pc has 
been found from fits to the position angles of polarize d starlight, giving a direction 
towards I, b = 38°, 23°, with uncertainties of ~ ±35° l|Frisch et aZ.ll201Clh . The ~ 70° 
angle between the magnetic field and upwind LSR CLIC direction is consistent with 
an ISMF compressed in an expanding superbubble shell. The ISMF direction indi- 
cated by the Interstellar Bo undary Explorer (IBEX) Ribbon arc center, £,b— 33° , 55° 
rtFrisch and McComasl[2010l . this volume) is consistent with the polarization direction 
to within uncertainties. A more distant measure of the ISMF is provided by fits to 
the Faraday rotation and dispersion measures for four pulsars, 150-300 pc away in the 
third galactic quadrant, which give a ISM F direction si milar to the polarization value, 
and indicate a field strength of ~ 3.3 fj,G (Sal vat i 201^). 

Not all nearby clouds are warm and diffuse. Towards the constellation of Leo, and 
within 12 pc, a tiny dense cold filamenta ry interstellar clou d with thickness < 0.4 pc 
has been identified in Na° absorption i Mever et qZjEooih . This cloud is similar to 
the tiny scale atomic structures (TSAS) that are observed throughout the ISM, with 
typical sizes 30 AU, volume densities ~ 10 cm , colu mn densities 10 — 10 cm~ , 
and thermal pressures P/k = nT = 10 4 - 10 6 cm" 3 K dStanimirovidl2009h . Once the 
solar motion is removed from cloud velocities (Figure 1, right), the small Blue cloud is 
also seen to form in a region where the LIC and G-cloud are colliding. Both the Leo 
cloud and the Blue cloud (e.g. HP 80007) coin cide with the ring of tiny dense cold 
clouds identified in H° 21-cm by dHaud l20inh. The CLIC is a decelerating flow of 
ISM, so that a scenario where the upwind G-cloud and downwind LIC are converging 
is consistent with the velocity data. Evidently tiny cold clouds can and have formed 
very close to the Sun through cloud collisions, so that the recent paleoheliosphere may 
well have included intervals where galactic cosmic ray modulation was minimal (see 
next section). 

Other possible types of nearby ISM include conductive, evaporative, boundaries 
between warm diffuse gas and hot plasma between the clouds. Unsaturated outflows 
of ISM from the clouds lead to sharp gradients in the ISM density, velocity, and tem- 
perature over small spatial scales. The outflow velocity can be up to 20 km s , with 
temperature variations of an order of magnitude, over distances < 0.5 pc (SF08). 



2.2 Configuration of Local ISM 

Our goal is to explore the connection between the paleoCISM and the paleomagnetic 
record on Earth. This connection is possible because of the heliosphere response to 
variations in the circumheliospheric ISM, and the cosmic ray flux at Earth is connected 
both to the heliosphere and neutrality of surrounding ISM. Optimally we need a model 
of the properties of nearby clouds, including cloud shapes, density, ionization, velocity, 
and homogeneity. Most of these data are unavailable, so we rely on simple assumptions 
for the cloud configuration and density. 

Discussions of the effect of the ISM on the heliosphere will be restricted to clouds 
that extend to within 5 pc of the Sun. For a typical Sun-cloud velocity of 17 km s , 
this would correspond to a look-back time of ~ 300, 000 years. The most complete data 
set available for determining cloud configurations are th e H° and D° data, which 
are free of uncertainties regarding elemental abundances (Rcdficld and Linsky||2004al ; 
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I Wood et ffldl2005h . Column densities for H° or D° can be transformed to distance 
scales characteristic of the neutral gas for some assumed volume density for H°. We 
will assume that nearby clouds have the same neutral density as the LIC, n(H°)=0.2 
cm~ 3 . 

Clouds near the Sun are typically identified through parsing absorption line veloc- 
ity data into different clouds by assuming rigid-body motion through space for each 
cloud. A numb er of studies provide data o n cloud velocities, but the most complete 
study is that of lRedfield and Linskvl l|2008l . RL08). However, one of the conclusions of 
this study, that the LIC does not surround the heliosphere, may change with better 
measurements of the velocity of interstellar He° inside of the heliosphere (Frisch and 
McComas, this volume). Distinct interstellar " clouds" have been found within 5 pc in 
a number of different studies (e.g. see review iFrisch et qilboilf ). We adopt the LIC 
and Blue cloud (B ) properties from RL08, and the G cloud (G) a nd Apex cloud (A ) 
fro m (IFris ch 2003); addit i onal co lumn density data is provided by I Wood et al\ (.2005) 
and lRedfield and Linskvl 1 200481 ) . The configuration of these clouds, within 60° of the 
galactic plane, is shown in Figure 1, right. The LSR velocity vectors for the clouds 
are shown, together with those of the CLIC and Sun. For the purposes of estimating 
cloud thicknesses from Ca + , the conversion factor of N( Ca + )/iV(HI) = 1.5 x 10 -8 
is adopted from a Aql and r\ UMa data. 

Estimates of the times it took (or will take) for the Sun to cross interstellar clouds 
observed towards several nearby stars are listed in Table 1, for the simplest assump- 
tion that the time is given by T — L/Vhc> f° r L — N(R°)/n(H ), heliocentric velocity 
VfjCi an d rl (H°)= 0.2 cm -3 . For stars near the upwind or downwind directions, this 
assumption should give realistic estimates as long as the cloud is homogeneous and 
neutrals and ions are well mixed. Obviously irregular cloud shapes, or cloud motions 
from unmeasured non-radial velocities, particularly for sightlines parallel to cloud sur- 
faces or for inhomogenous clouds, can affect these basic estimates. Also listed in Table 
1 are the angles between each star and the upwind directions of the vector CLIC, and 
vector LIC (from He°), velocities through the LSR. Whether clouds in the local ISM 
are filamentary or blobs affects the time inferred for the Sun to cross the cloud, partic- 
ularly for stars making a large angle with the upwind direction where radial velocities 
are small. 

The times in Table 1 are estimates only, since accurate dating of solar transitions of 
cloud boundaries requires data on the volume densities that are not generally available. 
If the clouds are in the rmal pressure equilibrium, t he cloud temperature variations of a 
factor of two or more (jRedfield and Linskvll2004bl ) suggest there are also variations in 
cloud densities, which can not be evaluated without information on cloud ionization. 
In addition, the structure of local ISM is resolved at very low spatial resolution, since 
current data only sample the sky at a spacing of ~ 1 star per 260 square-degrees (e.g. 
RL08). 



2.3 Possible Historical Changes in the Circumheliospheric ISM 

Stars located close to the downwind directions of cloud velocities will give the best 
estimates for entry times into the cloud, for simple geometries. For the LIC, downwind 
corresponds to small values for 0HeI- Estimates of the time when the Sun entered in 
the LIC will be based on interstellar data towards four stars, two near the downwind 
direction of the LIC (x 1 Ori and a Aur), and two at more oblique angles to the LIC 
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Table 1 Crossing times for Clouds Close to Sun C 1 ) 



Star 


I, 6,Dist,Cld 


Vkc> ^LSR 


N(H°) 


L, Cloud 


0CLICi "feel 


T, Crossing 










Thick. 


LSR, HC 


Time 




(deg,deg, pc) 


(km s- 1 ) 


(cm" 2 ) 


(pc) 


(deg) 


(years) 


a Cen 


316, -1, 1.3, G 


-18.0, -18.6 


17.61 


0.66 


20°, 130° 


(35,900) 


a CMa 


227, -9, 2.7, L 


19.6, 2.2 


17.60 


0.65 


106°, 42° 


32,200 


a CMa 


227, -9, 2.7, B 


13.7, -3.7 


17.40 


0.41 


106°, 42° 


29, 000 


a CMi 


214, 13, 3.5 


24.0, 10.4 


17.81 


1.86 


122°, 41° 


42,600 


a CMi 


214, 13, 3.5, L 


20.5, 6.5 


18.08 


1.38 


122°, 41° 


92,900 


70 Oph 


30, 12, 5.1, G 


-26.2, -9.3 


17.77 


0.95 


58°, 26° 


(670,000) 


70 Oph 


30, 12, 5.1, A 


-32.4, -15.5 


17.46 


0.47 


58°, 26° 


(530,000) 


a Aql 


48, -9, 5.1, G 


-18.1, -0.3 


17.91 


1.3 


72°, 50° 


(l.lefi) 


a Aql 


48, -9, 5.1, A 


-26.9, -10.7 


17.47 


0.48 


72°, 50° 


(670,000) 


X 1 Ori 


188, -3, 9, L 


22.3, 27.4 


17.93 


1.38 


146°, 13° 


60,500 


a Aur 


163, 4, 13, L 


21.8, 15.3 


18.24 


2.82 


171°, 28° 


126,000 


36 Oph 


358, 7, 6, G 


-28.2, -16.7 


17.85 


1.15 


26°, 170° 


(39,800) 



W Column 1 gives the star number in Figure 1, and star name. The galactic coordinates, 
distance, and cloud to which the component i s attributed , are listed in Column 2. G and 
A refer to the G-cloud and Apex cloud (fro m Frisch 200^), and L and B refer to the LIC 
and Blue clouds (from IRedfield and Linskvl [20081) . Colum n 3 gives the observ e d vel ocities 
of the interstellar absorption components in the solar (e.g., Rcdficld and Linskv (2008)) and 
LSR iner tial systems Co l umn 4 is the log of the cloud column density, from IWood et al\ 
(20051) or iHebrard et al\ | U999T) . Column 5 is the cloud thickness, AT(H°)/n(H°) calculated 
for volume density n(H°)= 0.2 cm -3 . Column 6 gives the angles #clic and 0Hel, which are 
the angles between the star and the LSR upwind CLIC direction (£, b = 335°, —7°), and the 
heliocentric (HC) downwind He° vector (£, b = 184°,— 15°), respectively. Column 7 gives a 
nominal crossing time for the cloud, calculated from the cloud thickness and Vhc- The times 
inside of parentheses are in the future. 



vector (a CMi, and a CMa). The Sun is assumed to be in the LIC. Because of the 
uncertainties in the ionization, volume density, magnetic pressure, an homogeneity 
of these clouds, the timescales below are quite uncertain. However they do show the 
plausibility of an interstellar effect on the cosmogenic isotope record. 

The star x 1 Ori, at 9 pc, is viewed through the heliosphere tail ((/>HcI~ L3°, Table 
1). The observed Vfjc = 22 km s _1 component corresponds to the LIC velocity. For 
the above premise that L =N(H°)/n(H°), the Sun would have entered the LIC ~ 
60, 500 year s ago. The relatively high turbulence of the LIC absorption lines towards 
X 1 Ori (e.g. IRedfield and Linskvf 2008") . however, suggests that two clouds with similar 
velocities could blend in velocity to form this component. If so, the Sun would have 
entered the LIC more recently. A more distant star in the downwind direction is a Aur 
(4>U c l= 28°). The higher LIC H° column densities suggest, instead, that the the Sun 
entered the ISM 126,000 years ago. The star a CMi is much closer to the Sun, but it 
is offset 41° from the downwind LIC direction (</>Hcl=41 ). The LIC component gives 
an entry time into the LIC of 92,900 years ago. 

The closest downwind star is Sirius (a CMa, 2.7 pc), with two clouds in front of 
it. Sirius is 106° from the downwind direction. The Sun would have entered the LIC 
component 32,200 years ago, and the Blue cloud ~ 29, 000 years before that. 

An alternate LIC encounter time follows from the assumption that the LIC is 
filamentary, with the filament oriented perpendicular to the LSR velocity and roughly 
parallel to the magnetic field direction. When distant diffuse ISM is spatially resolved, 
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filamentary structures are generally observed for neutral and ionized gas. Magnetic 
pressure creates filamentary ISM towards Loop I at distances of ~ 90 pc, so the evident 
relation betw een Loop I and nearby ISM suggest the LIC could also be filamentary. 
iFrischl (jl994T ) assumed a filamentary structure for the LIC, with the direction of the 
LSR cloud velocity perpendicular to the filament, and found the Sun may have entered 
the CISM gas anytime within the past 2,000-8,000 years, or even more recently. 

For the comparisons with the geologic radio isotope record in Section 4, and since 
the detailed structure of the downwind gas is not known, we adopt several time intervals 
as candidates for a solar transition between interstellar clouds. The primary downwind 
clouds of interest are the LIC, which is most likely presently surrounding the Sun, and 
the Blue cloud (Figure 1, right). These transitions are: (1) Sometime within the past 
8,000 years, as indicated by the a CMa sightline together with an assumed filamentary 
cloud structure. (2) Sometime within 10, 000 — 32, 000 years ago. (3) Before that the 
Sun was in the Blue cloud, with an entry time of ~ 61,000 years ago from the a 
CMa timescales. The Blue cloud should be more ionized than the LIC because it is 
less shielded from the primary source of local H-ionization, e CMa. This suggests that 
the Sun traveled through a cloud with higher ionization levels than the LIC, before 
it entered the LI C. This interpreta tion is also consistent with the local ionization 
gradient found by [ Wolff et all (|l999T ). (4) Any of the surfaces of these clouds may have 
a conductive interface. However the dumpiness of the local ISM, and the apparent 
deceleration evident in the CLIC velocities, suggest that the clumps of gas that form 
the CLIC have a cohesive origin and are close to each other. Transition times across 
0.5 pc cloud interfaces, where the relative Sun-cloud velocities may be larger by ~ 
20 km s , may be 10,000 years, or less, depending on the angle of the ISMF that 
inhibits conduction. (5) No tiny dense clouds have been directly detected nearby in the 
downwind direction, however the Blue cloud coincides with a ring of tiny cold dense 
clouds (Section I2.1|l . Typical sizes for these clouds, as small as 30 AU, indicate that 
such clouds would sweep past the heliosphere over timescales of a decade. 



3 The Heliosphere for Different Interstellar Environments 

The heliosphere results from a balance between solar wind and interstellar pressure. 
This balance sets the size of the heliosphere and determines the particle distributions 
and magnetic field configurations throughout, which in turn determine how galactic 
cosmic rays are being modulated during their passage through the heliosphere on their 
way to Earth. 

The inner boundary conditions of the heliosphere consist of the solar wind, im- 
printed with the 22 year magnetic activity cycle of the Sun, as a "point source" inflow. 
The time variability of the solar wind affects the large-scale heliosphere, including the 
heliosheath regions immediately inside and outside of the heliopause. The distance 
of the termination shock in the nose direction varies by about 10 AU with the solar 
cycle phase, and is clo ser during solar wind minima when ram pressures are lower 
(|Zank and Mullerll2003h . The reversals of the solar magnetic polarity every ~ 11 years 
between successive solar minima propagate to the inner heliosheath regions, and create 
bands of magnetic polarity t hat are swept around the flanks of the heliosphere with 
the subsonic solar wind fe.g.. |Pogorelov et all\2Q09l ). The outer boundary conditions of 
the heliosphere are dominated by the total pressure of the surrounding cloud, including 
magnetic pressure, the ram pressure of excluded ions, and a large fraction of the ram 
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Fig. 2 Left. The density distribution for the heliosphere when immersed in a moderately dense 
cloud, n(H°)~ 15 cm -3 , with a te mperature T=3,00 K, ion density of 0.2 cm -3 , and moving 
at the LIC velocity (Model 17 in IMiiller et qj]|2006f> . The heliopause is 34 AU from the Sun. 
Right. Spectra of galactic cosmic rays at 1 AU for three different interstellar clouds surrounding 
the heliosphere (ACRs are not included). LIC stands for the contemporary CISM; LB is the 
Local Bubble interior modeled as a 1.2 X 10 6 K fully ionized plasma with density 0.005 cm -3 ; 
and DC is a dense cloud with density 10 cm -3 , T=200 K, and relative ve locity of 25 km s -1 . 
Mode ls 1 and 2 utilize different cosmic ray modulation models (figure from lFlorinski and Zank 
12006ft ■ 



pressure of the neutrals which participate through charge-exchange with ions outside 
the heliopause; the pressure interior of the heliopause is modified by charge exchange 
as well. 

Configurations of the global heliosphere ha ve been modeled fo r a range of surround- 
ing interstellar cloud properties. In one study ( Miiller et fflZj20oih . the parameter space 



around the current interstellar cloud is sampled to establish empirical relations for lo- 
cations of heliopause, termination shock, and bow shock as a function of interstellar 
boundary values. This study confirms the overall connection of the size of the helio- 
sphere with the pressure balance argument. The interstellar ionization ratio weakly 
influences the downstream termination shock distance as well as the typical termina- 
tion shock strength. In relative terms, the termination s hock shape does not change in 
this particular region of parameter space. IMiiller et al. I (|2006h explore a wide range of 



possible interstellar environments, including cold dense clouds, hot tenuous, completely 
ionized clouds like that assumed for the Local Bubble, and systems on galactic paths 
that result in fast relative Sun-ISM velocities. In particular the relative velocity has 
a decisive influence on how much interstellar neutral material reaches the inner he- 
liosphere (filtration); the larger the relative velocity, the less filtration is occurring. If 
the Sun encounters a completely ionized ISM, such as the Local Bubble, GCR particle 
fluxes at Earth are lower than currently. With neutrals absent, none of the pressure 
balance modifications take place, and secondary particles like anomalous cosmic rays 
(ACR) do not exist. Lastly, if the heliosphere is embedded in a dense cloud, the he- 
liosphere is small in size, and particle fluxes at Earth rise substantially. Depending on 
the interstellar ram pressure, the heliosphere can easily get small enough for the Earth 
orbit to be partly in the inner heliosheath region. 
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lYeghikvan and Fahrl (|2006T l consider the passage of very dense molecular clouds 
over the heliosphere. The resulting heliospheres are very small, so that the orbit of 
Earth takes it through regions of interstellar gas outside the heliopause. The expected 
neutral fluxes at Earth are so high that changes for the ozone layer and other climate- 
related effects a re to be expected. Similarly, when a supernova shock front washes over 
the heliosphere ( Miiller et al\ I2OO9T ) . the heliosphere can be equally small (Figure [2^), 
with direct access of supernova material to the Earth atmosphere. In both cases of 
extremely sm all heliospheres, the ac cess of GCR to Earth is greatly enhanced. In a 
calculation bv lFlorinski et al\ l|2003T l. where a cold dense cloud results in a moderate 
size heliosphere of 23 AU, a two-fold flux of GCR, and a ten-fold flux of ACR is 
reported. Figure fjp illustrates this with calculated spectra for the Local Bubble case 
and for a dense cloud case, comparing them to the GCR fluxes of the contemporary 
heliosphere. 



4 Comparing past cloud transitions with the terrestrial radioisotope 
record 



Spallation of cosmic-rays in the terrestrial atmosphere creates the radioisotopes Be, 
36 CI, and 14 C that are used to date the repository geological archive. Both the geo- 
magnetic and he liosphere magnetic fields modulate the incident cosmic-ray flux (e.g. 
iMcC rackenll2004l . and the Beer and McCracken articles in this volume), including the 
ACR component. The well known anticorrelation between cosmic ray fluxes at 1 AU 
and the solar magnetic activity cycle (see articles by Leske and Mewaldt in this vol- 
ume) suggests that ISM-driven variations in the heliosphere may also be significant. 
Earlier interstellar explanations for the peaks in the concentration of 10 Be in the ice 
core record include the reduced modulation of co smic rays by a hel iosphere that has 
been compressed by a passing supernova shock ( Sonett et ql.lll987l and see Fig. 2), 



and higher GCR fluxes incident on the helios phere because the Sun traveled through 
a magnetic flux tube in the ISM (|Fnsch| j 



to climate cooling (e.g.. iKirkbv and Carslaw 



1997). High GCR fluxes have also been linked 



W- Hig h 
v 2006). 



The paleomagnetic field that is derived from geological data such as magnetic re- 
menance should, in principle, agree with the paleomagnetic field that has been derived 
from the radioisotope 10 Be and 36 C1 ice core data, and the 14 C tree-ring data. Recon- 
structing the geomagnetic field from cosmogenic isotope data requires an assumption 
about the flux and spec t ra of cosmic rays, and the production and dispersal of radioiso- 
topes. |Muscheler^FttI| ( 20051 . M05) compared the paleomagnetic field determined by 



these different techniques, and found several discrepancies that are not explained. Given 
the mystery of these anomalies, we postulate that temporal variations in the properties 
of the ISM shaping the heliosphere led to variations in the total cosmic-ray spectrum 
at the Earth. These variations have been overlooked as a factor in the geomagnetic 

time-line developed from cosmogenic isotopes, but need to taken into consid eration. 

Th e anomalies in the reconstructed paleomagnetic record mentioned in the lMuscheler et al\ 
(|2005l ) paper, and that are discussed here are: 1. The discrepancy in the timing of 
the maximum geomagnetic field dipole strengths (the "virtual axis dipole moment", 
VADM) during the late Holocene traced by the 10 Be versus 14 C records (see Fig. 6 in 
M05). 2. The interval of 18,000-34,000 years ago, where both 10 Be and 36 C1 VADM re- 
constructions show increased modulation compared to the remanence records, and the 
14 C data are poorly understood. 3. An anomaly 48,000 years ago where both 10 Be and 
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36 C1 showed increased modulation compared to the VADM record, and an anomaly 
58,000 years ago where only 36 CI shows increased modulation. 

Cosmic ray fluxes at Earth vary with the size and properties of the heliospheric 
modulation region. A second simple change in the CISM properties, that would neces- 
sarily affect cosmic ray fluxes at the Earth, would be spatial variations in the ionization 
levels of the interstellar gas. At 1 AU the cosmic-ray fluxes below ~ 200 MeV/nucleon 
are dominated by anomalous cosmic rays (ACRs). For instance, the ACR oxygen in- 
tensity of 10 MeV /nucleon is an order of magnitude larger than for GCR oxygen 
( Leske et qf]|2008h . ACRs form from neutral interstellar atoms that charge-exchange 
with the solar wind, creating pickup ions that are accelerated to cosmic ray energies. 
The ACRs will vanish if the surrounding ISM is ionized (with the exception of a very 
minor contribution arising from dust-generated pickup ions). These spatial variations 
translate to temporal variations in the 1 AU cosmic ray fluxes because of the relative 
Sun-cloud motion. 

Models of radioisotope production include the GCR sources, and cosmic rays from 
energetic solar flares. However, they do not distinguish ACRs as a source with a distinct 
origin and spectrum, and one that is likely to vary as the Sun encounters different 
clouds. Radioisoto pe production has b een calculated for solar cosmic rays at similar 
energies as ACRs \ Webber et qZj|2007h . bu t not for the ACR s pectrum itself. ACRs 
are trapped in the radiation belts of Earth i Mewaldt et aZ.1119981 ). 

The primary question then becomes: How would the 14 C, 10 Be, and 36 C1 production 
change as the ACR component vanishes? Three possible effects may contribute. The 
first is that 14 C is formed by thermalized neutrons near the top of the atmosphere, 
and th e long storage time (> 0.1 year) of ACRs trapped in the radiation belts of 
Earth I Mewaldt et a/.lfl998h may multiply the production of 14 C because of the long 
exposure times of nitrogen and carbon compounds in the atmosphere compared to the 
direct production of cosmogenic isotopes formed by spallation. An opposite effect of 
ACRs on cosmogenic isotope production is suggested by the yield functions of 10 Be and 
14 C as a fu nction of the e nergy of the incident cosmic ray proton, which are shown in 
Figure 6 in lUsoskinl (|20Q8h . The production of 10 Be is relatively more efficient at ACR 
energies then the production of 14 C, based on comparisons of the ratios of production 
at ~ 200 MeV (ACR energies) and ~ 20 GeV (GCR energies). These ratios for 10 Be 
and 14 C are ~ 0.03 and ~ 0.001, respectively. These yield functions would then suggest 
that the production of 14 C is less sensitive to the ACR component of cosmic-ray fluxes 
than is the production of 10 Be. A third wild-card possibility is that rapid variations of 
the geomagnetic field may reduce coupling between the ACRs and the radiation belts, 
so that ACRs have the same access to the atmosphere as GCRs, reducing any effect 
of storage of ACRs in the radiation belts. In the absence of a detailed understanding 
of the effect of ACRs on cosmogenic isotope production, the discussions below linking 
cloud transitions to the geological radioisotope data are highly speculative. 



4.1 Late Holocene Discrepancies and Anomalous Cosmic Rays 

The VADM paleomagnetic field reconstructed from 10 Be records peaks at 2,000 years 
BP, which is one millennium after the peak determined from 14 C at 3,000 years BP (Fig. 
6 in M05). At the CHISM velocity, 3,000 years corresponds to a distance of 0.08 pc, 
which is a plausible value for the distance to the cloud surface in the a CMa direction. 
The 2,000-3,000 time interval is within the uncertainties of the solar entry into the LIC 
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as inferred from the a CMa sightline bv lFrischllj 19941 ). based on a filamentary structure 
for the L1C. A LIC filamentary structure would be consistent with Figure 1, right, only 
if the density structure of the ISM moving at the LIC velocity is inhomogeneous. 

Since VADM variations anticorrelate with radioisotope fluxes, th e 14 C fluxes de- 
crease d a millennium before the 10 Be fluxes decreased in the Holocene. lMuscheler et all 
( 20051 ) postulated that this discrepancy was due to the carbon cycle, because the 10 Be 
and archeomagnetic field determinations generally agree. 

If the discrepancy is due, instead, to variations in the ISM, we suggest that ACRs 
may be the culprit. The three possible influences of ACRs on the cosmogenic isotope 
record leads to three different explanations for the anomaly. Prior to entering the CISM, 
the Sun may have traversed a fully ionized cloud interface with no neutrals, no pickup 
ions or ACRs, with a larger relative Sun-cloud velocity leading to higher interstellar 
ram pressures (see Fig. 2 in SF08), and with reduced mass- loading of the solar wind by 
pickup ions. (1) If the enhanced yield of 10 Be production at low energies is important, 
the absence of ACRs at 1 AU would increase 14 C relative to 10 Be, which is not seen. 
(2) An opposite result follows if the long exposure times for 14 C formation from ACRs 
in the radiation belts is important, in which case an ionized ambient ISM would prefer- 
entially reduce 14 C levels because of the concordant lack of 14 C enhancement from the 
multiplying effect of the storage of ACRs in the magnetosphere. Once the Sun enters 
the main part of the cloud, which is ~ 75% neutral, the low energy ACR component 
is restored and the balance between the synchronization of the 10 Be, 36 C1, and 14 C 
formation was restored along with the ACRs. (3) For the third alternative, if disrup- 
tion of the terrestial magnetic field affects the 1 C production only, 10 Be would stay 
coupled to the geomagnetic field while the timescale for 14 C dating would be disrupted. 
The tight coupling between 10 Be and the geomagnetic field suggests that the overall 
heliospheric modulation of the GCR fluxes did not change, so that the discrepancy 
could be due to the effect that the disrupted geomagnetic field has on the flux of ACRs 
into the atmosphere. 



4.2 Cloud Crossings 

During the period 18,000-34,000 years ago, both the 10 Be and 36 C1 VADM recon- 
structions deviate from the archeomagnetic remanence records. The increased VADM 
reconstructed from 10 Be and 36 C1 suggests reduced fluxes below t hat explained by 
the ar cheomagnetic record, or extra heliosphere modulation of GCRs. iMuscheler et al\ 
(2005) suggested that possible changes in levels of solar magnetic activity may be re- 
sponsible for this difference. We speculate that this interval instead corresponds to an 
encounter between the Sun and slower clumpy ISM beyond the CISM in the downwind 
direction that forms the second clouds observe d towards a CMa (the Blue cloud, Figure 
1 right) and towards a CMi (the Aur cloud in iRedfield and Linskvll2008l ). both within 
3.5 pc. Cloud velocities alone affect the GCR modulation because the heliosphere di- 
mensions increase as the ram pressure of the ISM decreases. The decreased 10 Be and 
36 C1 fluxes in this interval would then suggest a larger heliosphere modulation region, 
that still contains a solar wind mass-loaded with pickup ions from interstellar neutrals. 
Cloud-Sun relative velocities of 26 km s _1 today, versus 21-22 km s _1 for x 1 Ori and 
a Aur, and 14 km s _1 for the Blue cloud towards a CMa (Table 1), suggest that cloud 
ram pressure could have been a factor of ~ 2 — 3 lower in the past than today, giving 
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a larger heliosphere and larger GCR modulation region that are consistent with the 
decreased fluxes found by M05. 

Additional discrepancies, about 46,000-48,000 and 58,000 years ago, are seen be- 
tween the radio isotope paleomagnetic and remanence paleomagnetic records. If the 
Sun entered the neutral portion of the LIC ~ 60, 500 years ago, as indicated by x 
Ori located in the downwind direction (Table 1), the VADM discrepancy suggests that 
increased cosmic ray modulation in the LIC, or decreased ACR production, could have 
produced the discrepancy. A slightly higher density, n(H°)~ 0.25 cm~ 3 , would match 
the solar entry into the LIC to 48,000 years ago discrepancy. 

5 Future: Improving Comparisons between Local ISM Structure and 
Paleomagnetic Records 

Comparisons between the structure and kinematics of interstellar clouds near the Sun 
and the paleomagnetic records deduced from radio isotopes such as 10 Be, 36 CI, and 14 C, 
are based on scanty knowledge of details of the three-dimensional spatial distribution 
and configuration of nearby interstellar clouds. Volume densities of both neutrals and 
electrons in the nearest interstellar clouds are required to improve maps of the ISM 
distribution. Some insight could be gleaned from photoionization models, as done for 
the CISM; however high-sensitivity ultraviolet data are required for realistic models. 

The most significant open question is whether ACR variations arc sufficient to 
require consideration as a separate component that contributes to the production of 
the cosmogenic isotopes, and that varies separately from the GCR component as the 
heliosphere is modulated by the ISM. One aspect of answering such questions will be to 
expand calculations of the fluxes of galactic and anomalous cosmic rays at the Earth 
beyond the limited number of heliosphere configurations now available. More study 
of the roles of ACRs versus higher GCRs as source populations of the radio isotopes 
are needed, since differences between the 14 C, 10 Be, 36 C1 records are mandated if 
ACRs are a factor in 14 C production rates. A better understanding of the physical 
properties and configuration of ISM close to the Sun is also required for accurate 
comparisons between the timelines of the cosmogenic isotope and heliosphere boundary 
conditions. Nevertheless, we conclude that in principle the attenuation of GCRs by 
an ISM-modulated heliosphere, and perhaps the presence or absence of ACRs, are 
capable of accounting for differences between the paleomagnetic record determined 
from cosmogenic isotopes versus remnance data. 
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